Comment on "New modes of halo excitations in the 6 He nucleus" 
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We try to explain the differences in the He dipole strength 
function in jl| and Q| . We perform the full basis calculation of 
the strength function with the same renormalized interaction 
as in Q and show that the size of the basis, needed for con- 
verged calculations of the 6 He continuum spectrum, is much 
larger than that for the discrete spectrum. The renormalized 
interaction of therefore cannot be used for the continuum 
spectrum calculations with the same basis as for the ground 
state. 



In a recent paper Q the continuum properties of 
6 He were discussed within a three body nna model. 
The strength functions were computed using the hyper- 
spherical harmonics approach. The authors argued that 
the main features of the continuum spectrum of the 
nna system can be obtained using only a few hyper- 
spherical harmonics with the hyper-spherical quantum 
number K < 6, if one simultaneously renormalizes the 
na interaction by introducing additional attraction. 

However the dipole strength function in |l]] differs sig- 
nificantly from our recent calculation where we used 
the same three-body model. The purpose of this com- 
ment is to explain these differences. 

We show that the continuum state calculations of the 
nna system demand a much larger basis, K < 100, com- 
pared to the bound state. In such case the simple pre- 
scription to renormalize the na interaction aiming at the 
bound state might not be adequate for the continuum. 

The method we use is, in principle, similar to that 
of - the angular part of the wave function is ex- 
panded in terms of the hyper-spherical harmonics and 
then the resulting system of coupled hyper-radial equa- 
tions is solved. However our approach has the follow- 
ing modifications ||: i) We use Faddeev rather than the 
Schrodinger equation; ii) Prior to solving the hyper-radial 
equations we calculate the eigenvalues A„, as functions 
of the hyper-radius p, of the angular part of the Faddeev 
operator. The quantities A„(p)p~ 2 then serve as effec- 
tive potentials in the radial equations; iii) We use an- 
alytic large- distance solutions of the Faddeev equations 
for short range potentials Q which greatly enhances the 
method and allows fully converged calculations. 

The na interaction, originally fit to the scattering 
data, was somewhat modified in JjLjsB in order to com- 
pensate for the limited basis, 

Vna{ r ) — exp (— r 2 /2.35 2 ) x 
(50 P s - 48.3 P p - 23.0 P d - 11.71 1 • s). 
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FIG. 1. The lowest angular eigenvalue for 6 He (J" = 1~) 
calculated with different K max . Also shown is the eigenvalue 
obtained from the asymptotic large distance expansion of the 
Faddeev equations §4|. 

All lengths are in fm and energies in MeV. Due to 
the over attraction this interaction fails to reproduce the 
experimental na phase shifts. 

For the nn interaction we use a simple Gaussian [Q, 
similar to the the one used in |1. The potential repro- 
duces the experimental effective range in s-wave and scat- 
tering lengths in s- and p- waves, 



V nn (r) =exp(-r 2 /1.8 2 ) x 

(2.92 + 45.22 Sis 2 - 12.08 1 • (si + s 2 ) 



(2) 



26.85 S 12 ), 



where si, s 2 are the spins of the neutrons, S12 is the usual 
tensor operator. 

In Fig. 1 we illustrate the convergence of the lowest an- 
gular eigenvalue using, as examples, calculations with 



K. 



9, K„ 



Here Pi is a projection operator on the corresponding 
orbital momentum state 1 and s is the neutron spin. 



19, and K max = 99. Obviously the 
basis limited to K < 9, used in Q, is by far not con- 
verged. It is only the curve for K max sa 100 that finally 
converged towards the large-distance asymptotics. This 
long tail of the eigenvalue arises, in the present case, due 
to large scattering length in the nn subsystem [Q. The 
scale of convergence towards the asymptotic behavior is 
then around 20 fm. 

The convergence problems for the lowest effective po- 
tential is illustrated in Fig. 2. The converged curve has 
the minimum value of about 0.15 MeV and the low (0.3 
MeV) barrier extending into the asymptotic region begin- 
ning at about 25 fm. In contrast the curve for K max = 9 
has minimum and maximum values at 0.7 MeV and 1.4 
MeV, respectively. A more realistic potential calculated 
with a correct na interaction , is accidentally close to 



the curve for K n 
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FIG. 2. The lowest effective potential (Ai (p) + 15/4)p -2 x 
% 2 /(2m) for 6 He (J™ = 1~) calculated with different K max 
(m is the nucleon mass). 

Note that that the K < 9 basis seems to be sufficient 
for the + ground state with the renormalized interaction 
jlj . This shows that the renormalization technique needs 
different bases for the ground state and the continuum 
state calculations. 

The dipole strength function was calculated in 
with a basis size corresponding to K < 6. In Fig. 3 we 
show the results of our calculations of the dipole strength 
function for different K max . Our K max — 9 curve agrees 
quite well with pL Increasing K max causes the main 
peak to shift towards lower energies with the second peak 
forming and also moving towards lower energies. The re- 
sult for K max = 9 is rather smooth and deviates rela- 
tively little from the corresponding result obtained with 
plane waves for the 1 _ continuum wave function. This 
leads to a presumably misleading impression that there 
is little structure in the 1~ continuum. The other two 
curves, K max =19 and K max =99, on the contrary show 
visible concentration of the strength at small energies. 
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FIG. 3. The dipole strength function S(E) = ^ n x 
|(n_E||M(i51)||(;s)| 2 as a function of the three-body energy E 
for the transition from the ground state (0 + ) to the contin- 
uum state (1~) of 6 He calculated with different K ma x- The 
curve from [jlj] is given in arbitrary units. 



We have to note, however, that due to the additional 
attraction introduced in the na interaction in ]jj the con- 
verged dipole strength function here is excessively shifted 
towards lower energies and therefore is not realistic. The 
realistic dipole strength function M is close to the curve 
Kmax = 19. In this case the limited basis is approx- 
imately compensated by the over-attractive potential. 
Note that for the + bound state calculations the cor- 
responding compensatory basis was K max = 6 |Q. This 
shows that the continuum spectrum demands more care- 
ful and accurate calculations. 

In conclusion, in order to explain the differences 
between the dipole strength functions from |Q and Q 
we have performed a fully converged calculation of the 
= 1" continuum spectrum of 6 He using the same 
renormalized na interaction as in Q. We have shown 
that the number of basis functions needed for converged 
calculations of the 6 He continuum spectrum is much 
larger than that for the discrete spectrum. We conclude 
that the technique of jij to renormalize the na interaction 
to compensate for the limited basis becomes ambiguous 
when applied to the continuum spectrum. The reason for 
the differences is therefore that this renormalized inter- 
action is not appropriate for continuum calculations with 
the same basis size as for the ground state. 
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